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ABSTRACT 

We present a Spitzer InfraRed Spectrometer (IRS) search for 10-36 /im molecular emission from a 
large sample of protoplanetary disks, including lines from H2O, OH, C2H2, HCN and CO2. This paper 
describes the sample and data processing and derives the detection rate of mid-infrared molecular 
emission as a function of stellar mass. The sample covers a range of spectral type from early M to A, 
and is supplemented by archival spectra of disks around A and B stars. It is drawn from a variety of 
nearby star forming regions, including Ophiuchus, Lupus and Chamaeleon. Spectra showing strong 
emission lines are used to identify which lines are the best tracers of various physical and chemical 
conditions within the disks. In total, we identify 22 T Tauri stars with strong mid-infrared II2O 
emission. Integrated water line luminosities, where water vapor is detected, range from 5 x 10^'* to 
9 X 10^"^ Lq, likely making water the dominant line coolant of inner disk surfaces in classical T Tauri 
stars. None of the 5 transitional disks in the sample show detectable gaseous molecular emission with 
Spitzer upper limits at the 1% level in terms of line-to-continuum ratios (apart from II2), but the 
sample is too small to conclude whether this is a general property of transitional disks. We find a 
strong dependence on detection rate with spectral type; no disks around our sample of 25 A and B 
stars were found to exhibit water emission, down to 1-2% line-to-continuum ratios, in the mid-infrared, 
while almost 2/3 of the disks around K stars show sufficiently intense water emission to be detected 
by Spitzer. Some Herbig Ae/Be stars show tentative II2O/OII emission features beyond 20 /im at the 
1-2% level, however, and one of them shows CO2 in emission. We argue that the observed differences 
between T Tauri disks and Herbig Ae/Be disks is due to a difference in excitation and/or chemistry 
depending on spectral type and suggest that photochemistry may be playing an important role in the 
observable characteristics of mid-infrared molecular line emission from protoplanetary disks. 

Subject headings: astrochemistry - planetary systems: protoplanetary disks - stars: pre-main sequence 



1. INTRODUCTION 

Molecular volatiles, including H2O, CO2 and CH4, 
among others, are thought to play a central role in the 
early evolution of planetary systems. Evidence collected 
from Solar System material has shown that the Solar 
Nebula was characterized by a diverse chemistry during 
the process of planet formation. This chemistry played 
a pivotal role in the formation of giant planets by pro- 
viding a reservoir of volatile solids to aid in a rapid 
build-up of planetary cores. It has long been thought 
that the volatiles that formed in the nebula seeded the 
surfaces and atmospheres of the terrestrial planets af- 
ter their formation, the so-called "late veneer" model. 
The mechanism for delivering a veneer of water to an 
otherwise dry Earth is controversial, although primitive 
carbonaceous chondrites provide the best match to the 
D/H ratio of the Ear th's oceans (for an overview, see 
iMorbidelli et al.l [20001 ). Alternatively, it has been sug- 
gested that Earth's water was present at earlier stages, 
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having been chemisorbed to small dust grains at temper- 
atures too high for the formation of ices — a scenario pos- 
sible if the Earth-forming olig archs formed in a disk env i- 
ronment rich in water vapor (jMurahdharan et al.l l2008'). 

Because these models are based on evidence gathered 
from ancient Solar System material, very little is known 
about how volatiles evolve and influence planet forma- 
tion in current analogs to the Solar Nebula — the so- 
called protoplanetary disks. The material that formed 
the present-day Solar System originated in feeding zones 
located at distances within ^20 AU from the Sun, known 
putatively as the •planet-forming region. It is a natural 
conjecture that solar nebula analogs can be found among 
the zoo of protoplanetary disks around young stars. 

Theory predicts that water plays a role in the physics of 
planet formation (in addition to the chemistry) . Because 
0-, C- and N-based volatiles, with H2O likely being the 
most abundant, account for half the condensible mass 
in a protoplanetary disk, the formation, transport and 
phase changes of these molecular species can profoundly 
affect the evolution of protoplanetary disks. First, theory 
predicts that the dynamic interplay between, in particu- 
lar, water vapor and the buildup of planetesimals is cru- 
cial for planet formation. Processes such as freeze-out, 
gas diffusion and inwards gas-drag migration of solids 
may act to concentrate large amounts of ice and water 
vapor close t o the midplane snow line, facilitating planet 
growth fe.g.. iStevenson fc Luninelll988t [Johansen et al.l 
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[2007t iDodson-Robinson eral|[200l iCieslal [20091 ) . It is 
still an open question how the inward and outward wa- 
ter transport processes balance, but it seems likely that 
strong abundance gradients may be present in the inner 
10 AU of protoplanetary disks, in both the radial and 
vertical directions. Such molecular abundance structures 
— not only for water, but for most molecular species — 
should be observable. 

The planet-forming region represents a chemical en- 
vironment very different than those studied in the in- 
terstellar medium, molecular clouds and even the outer 
regions (>100 AU) of the same disks. The very high den- 
sities, temperatures, radiation fields and, importantly, 
short dynamical and chemical time scales in such re- 
gions generate conditions that may be more easily com- 
pared to the chem istry of planetary atmospheres (e.g., 
iWoitke et aDl2009( ). In direct comparison, models for 
planetary atmospheres show that complexities are high 
and predictabilities low, making the acquisition of empir- 
ical data crucial for understanding inner disk chemistry 
and disk "climate" . 

The recent detections of a large number of rotational 
emission lines due to warm water vapor, in addition to 
the rovibrational Q-branchcs of HCN, C2H2 and CO2 
in three protoplanetary disks — AA Tau, AS2 05N and 
DR Tau (|Carr fc Naiital[2008t ISalvk et al.ll200l — have 
prompted a number of questions: How common is this 
emission? If a large number of protoplanetary disks ex- 
hibit strong molecular emission in the mid-infrared, the 
observed lines have the potential to be a unique tracer 
of inner disk chemistry as well as disk evolution and 
planet formation processes. Are differences from disk 
to disk present? If many disks exhibit similar molecu- 
lar signatures, they will serve as good tracers of physi- 
cal and chemical conditions or disk evolution. A funda- 
mental question is therefore which physical and chemical 
processes result in the presence or absence of molecular 
emission tracers. If lines are absent in some disks, for ex- 
ample, is it due to abundance or excitation differences? 

We have collected a database of high quality Spitzer 
spectra aimed at the analysis of gas-phase emission lines, 
based on both new, dedicated observations, as well as 
on re-processed archival data. Given the large amount 
of information contained in each spectrum, we anticipate 
presenting a series of papers focusing on different aspects 
of the data. In this first manuscript, we present the data 
and use a sample of ~75 protoplanetary disks spanning 
a wide range in stellar effective temperature to demon- 
strate that the presence of complex molecular emission, 
including that of water vapor, is indeed a common prop- 
erty of protoplanetary disks around stars of solar mass 
and lighter. We also show that disks around A and B 
stars do not show strong molecular emission, as defined 
by having line-to-continuum ratios lower than their T 
Tauri counterparts by an order of magnitude or more. 
A companion paper presents a detailed analysis of the 
molecular emission spectra, including a comparison to 
emission models and a detailed study of the full range of 
molecular features, including HCN, C2H2, CO2 and OH 
features (Salyk et al. 2010, herafter Paper II). Future 
work will describe a companion study of T Tauri disks of 
the Taurus star forming region (Carr et al., Najita et al., 
in prep), as well as more detailed modeling and compar- 
ative studies of stellar properties such as UV and X-ray 



fluxes, accretion rates, etc. 

2. OBSERVATIONS 

2.1. Sample selection 

The young star-fdisk systems examined in this pa- 
per are drawn from several Spitzer-IRS short-high (SH) 
and long-high (LH) observing programs. Most of the 
Chamaeleon and Lupus sources were observed as part 
of a deep T Tauri star survey (PID 50641, PI. J. Carr, 
hereinafter GO-5) that was designed to obtain very high 
S/N spectra by maximizing redundancy and using dedi- 
cated background observations for each target. The sub- 
set of the GO-5 sample presented here was constructed 
by selecting disks spanning the known range in disk in- 
frared colors and stellar X-ray luminosity, in particular 
those previously o bserved in the the cores to disks (c2d) 
Legacy program ()Evans et al.l 120031 ) that showed some 
hint of water emission (but which were not observed with 
sufficient redundancy and high enough S/N to produce 
firm detections) . Another subset of targets from the high 
quality survey focused on the Taurus cluster, and will be 
the subject of a separate study (Carr et al., in prep). 
The sample covered here is summarized in Tables [T] and 

m 

The non- Taurus 'high-quality' sample is supplemented 
by additional archival c2d spectra; and for comparison 
with the T Tauri stars that dominate the c2d selection, 
spectra from t he Herbig Ae/ Bc star survey by J. Bouw- 
man (PID 3470 lBoersma et al. 2008) were extracted from 
the Spitzer archive and reduced using the same proce- 
dure as the remaining sample. This ensures that the 
sample includes a significant number of disks associ- 
ated with stars of spectral types spanning from M to 
B. The full sample includes disks from the young clus- 
ters in Perseus, Taurus, Chamaeleon, Lupus, Ophiuchus 
and Serpens, but is not complete in any strict sense. 
For instance, the sample represents a selection of the 
"best" available Spitzer spectra, and thus tends to be 
biased toward brighter, isolated, disks with low back- 
gr ound emission (the c2d selection criteria are described 
in lEvans et al.ll2003l ). Nevertheless, given the size of the 
sample (75 disks), relative to the total number of class 
II disks brighter than ~100 mJy at 8/im in the major 
nearby star forming cloud s surveyed by the c2d program 
(~ 120. lEvans et alJl2009( ). there is little room for strong 
biases, except for those imposed by the IRS sensitivity 
limit. A disk with a fiat spectrum and a flux of 100 mJy 
at 8 fj,Ta roughly corresponds to a star with mass ~0.1 
Mq at the distance of Ophiuchus (125 pc), and 0.5 Mp, 
at the distance of Serpens (41 5 pclDz ib et al.''2010'). us- 
ing the evolutionary tracks of iSiess et al. (2000). 

Figure [1] shows the distribution of spectral in- 
dices ^1 ri jlm,—^^ jim. of the s ample, a>s d efined in 
iKessler-Silacci et al.l H2006f ) and iFurlan et al.l (|2006[ ). A 
few sources are labeled, for reference, including the tran- 
sitional disks T Cha, LkHa 330, HD 135344B and SR 21. 
It can be seen that the sample spans the range from disks 
with slopes that are nearly photospheric, fii3^m_3ipm'^- 
2, to transitional disks with large inner holes and rising 
fluxes into the far-IR ni3^,n-3ipm'^2. The majority of 
the disks center around 71.13 ^m_3i^m'-^-l ^ 0, a distribu- 
tion quite similar to that of the general disk populations 
of the included star forming regions (Furlan et al...2009.) . 
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Fig. 1. — Distribution of mid-infrared spectral indices and 31 
/im fiuxes of the sample. The stars indicate T Tauri stars (spectral 
type <F), while the squares denote Herbig Ae/Be stars. 

The Herbig Ae/Be stars tend to be brighter than the T 
Tauri stars, but not enough to reflect their much higher 
luminosities - an indication that the Herbig Ae/Be stars 
arc, on average, located at larger distances than the T 
Tauri stars. The sample also shows a deficit of tran- 
sitional disks with low 31 /im flux. This may, to some 
extent, reflect an intrinsic property of transitional disks, 
although one canonical transitional disk that is not in 
our sampl e — GM Aur — o nly has a 31 /im brightness 
of - 1 Jy ()Furlan et al.l[2006[ ). 

2.2. Data reduction 

The Spitzer-IRS spectra are reduced using IDL scripts 
optimized for deep integrations with high redundancy, 
i.e., up to 56 individual spectral frames for our data set, 
and dedicated backgroun d observations. The p rocedure 
is similar to that used by iCarr fc Najital (|2008[ ) for AA 
Tau, and the resulting spectra have been benchmarked 
to those results. The reduction begins with the droop 
frames. On-source spectra as well as background expo- 
sures are co-added. To avoid inappropriately weighting 
pixels at the edges of the entrance slit, the flat fleld is di- 
vided by low-order polynomial fits in both the dispersion 
and cross-dispersion directions. This produces a flat fleld 
that corrects the pixel-to-pixel response only. The back- 
ground exposure is used to detect rogue pixels by flagging 
pixels with values that are more than 2a from the mean 
of that pixel. Next, the background frame is subtracted 
from the on-source frames (two nod positions) and the 
bad pixels are linearly interpolated in the dispersion di- 
rection. The ability to detect transient bad pixels by 
using a highly redundant (i.e. including many individ- 
ual readouts) background observation taken essentially 
at the same time as the pointed observation produces 
excellent results. An example LH frame, before and af- 
ter bad pixel cleaning, is shown in Figure [5J 

Spectra are extracted from the 2D co-adde d spectral 
images using optimal extraction (|Hornelfl986f ). A selec- 
tion of 15 standard star observations of S Dra and ^ Dra 
were retrieved from the archive and reduced in the same 
manner, with the same flat fleld and extraction aper- 
tures, essentially producing a data base of spectral re- 
sponse functions (SRFs). Because the spectral response 
function depends on the location of the target in the 
slit, a data base of SRFs allows a selection of those that 
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Fig. 2. — Quality of the Long-High (LH) cleaning procedure and 
background subtraction. Dedicated background frames as well as 
high frame redundancy are especially critical for high SNR LH 
observations, by enabling an efficient removal of bad pixels. 

best match a given science observation. The best 6-12 
standard stars were chosen by minimizing the noise at 
10-11 /im (SH orders 19 and 20) and a flux calibrated 
spectrum is produced for each standard star observation. 
The flux calibrated spectra arc carefully defringed using 
IRSFRINGE (Lahuis et al. 2007) ensuring that the real 
structure from the complex water emission spectrum is 
not affected. Setting the tolerance of a defringer too low 
will essentially result in the data being processed by a 
low-pass fllter, removing high frequency structure, such 
as that produced by the densely packed water spectrum. 
The behavior of IRSFRINGE was tested by visually com- 
paring the flnal product with spectra that were not de- 
fringed. In practice, most orders required no defringing 
at all, with most of the fringing pattern being removed 
by the SRFs. The orders are then combined using a 
weighted average with the blaze function as weights. No 
additional relative scaling of the orders is necessary, as 
the match is generally within the noise. Finally, the set 
of standard star divided spectra that contain the least 
amount of residual fringes are co-added to produce a fl- 
nal spectrum. 

For the set of spectra that do not have dedicated back- 
ground observations, such as the suite of c2d spectra, a 
background observation taken as close in time as possi- 
ble (<6 months, typically 1-3 months) is used instead. 
Because no dedicated background observations were ob- 
tained for the c2d spectra, a few sources suffer from an 
incomplete background subtraction. However, given the 
spatial undersampling of Spitzer high resolution spectra, 
it is difficult to subtract the background without intro- 
ducing strong artifacts. 

3. THE MID-INFRARED MOLECULAR SPECTRUM OF A 
TYPICAL PROTOPLANETARY DISK 

The molecular emission spectrum from an optically 
thick protoplanetary disk, where detected, is charac- 
terized by >100 line complexes present throughout the 
Spitzer high resolution 10 - 36 /im wavelength range. 
The vast majority of lines are due to rotational transi- 
tions of the main water isotopologue, H2®0. Recall that 
this asymmetric top molecule gives rise to three rota- 
tional quantum numb ers. We use the HITRAN catalog 
(|Rothman et al.l |2005[ ) for line identification and adopt 
the standard notation, such that a rotational level is de- 
fined by JKaKc- In general, the level energy increases, 
within a J state, with the difference Ka — Kc, with a few 
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Not e. — REFERENCES - (1) Gucnthcr ct al. (2007); (2) Lcvrcault (1988); (3) Fernandez ct al. (1995); (4) Casali fc Eiroal 
([1993): (5) [ Kenyon fc Hartmann (1995); (6) Nguyen et al, t2009): (7) Hartmann et al. (1998) (8) Aleala et al. (1993); (9) 
Sehisano et a l. (2009); (10) Hughes & Hartigan (19 92); (1 1) iCoulson fc Walther (1995); (12) Stempels et al. (2007); (13) 
Stempels fc >iskunoy (2003); fl4) Seperuelo Duarte et al.l 1I2OO 8): (15) Hughes c t al.. (.199 4) ; (16 ) C_oh cn fc Kuhi (1979); 
(17) Martin (1998); (18) Natt a et^ ( I2OO6I); (19) Wilk ing et al.l (120051): (20) IPrato" et "Sfl (|2003) ; (2 iVlLee 6^0] (119941 ): 
(22) Bouvier fc Appenzcllcr (19921); (23) iRvdgrcn (19801): (24) IGrankin et al.l ( 120071) ; (25T Trorres et al.l ( 120061) ; (26) 
iStephenson fc Sanduleah 019771) 

Type of background observation: 1 - dedicated background observation obtained at the same time at the on-source observation. 
2 - "best effort" archival background observation. 
^ From lCarr"07^aiita (,2008 ) 
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TABLE 2 

Sample of protoplanetary disks around Herbig Ae/Be stars (>A) 



Source name Sp. Type Distance SH int. time LH int. time BG obs.^ AOR Obs. program 
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A3-A5i'8 


200 


4 


X 6.3 


4 X 6.3 


1 


rll001088 


3470 


HD 244604 


AOi'S 


400 


6 


X 6.3 


8 X 14.7 


1 


rll001344 


3470 


HD 36917 


B9.5/A0.52'* 


400 


4 


X 6.3 


6 X 14.7 


1 


rll001600 


3470 


HD 37258 


Al-A2^'® 


506 


8 


X 6.3 


4 X 61 


1 


rl0998784 


3470 


BF Ori 


A54 


400 


2 


X 31.5 


2 X 61 


2 


r5638144 


c2d 


HD 37357 


A05.8 


506 


8 


X 6.3 


6 X 14.7 


1 


rll001856 


3470 


HD 37411 


B95,8 


506 


6 


X 31.5 


8 X 14.7 


1 


rll002112 


3470 


RR Tau 


A0« 


2000 


2 


X 31.5 


2 X 61 


2 


r5638400 


c2d 


HD 37806 


B9-A2i'8 


473 


4 


X 6.3 


6 X 6.3 


1 


rll002368 


3470 


HD 38087 


B58 


473 


4 


X 31.5 


8 X 14.7 


1 


rll002624 


3470 


HD 38120 


B98 


506 


4 


X 6.3 


4 X 6.3 


1 


rll002880 


3470 


HD 50138 


B88 


290 


4 


X 6.3 


4 X 6.3 


1 


rl 1003648 


3470 


HD 72106 


AO* 


290 


8 


X 6.3 


6 X 14.7 


1 


rll004416 


3470 


HD 95881 


A18 


118 


4 


X 6.3 


8 X 6.3 


1 


rl 1004928 


3470 


HD 98922 


B98 


1000 


2 


X 6.3 


2 X 14.7 


2 


r5640704 


c2d 


HD 101412 


AO* 


160 


2 


X 31.5 


2 X 61 


2 


r5640960 


c2d 


HD 144668 


A7» 


208 


4 


X 6.3 


4 X 6.3 


1 


rll005952 


3470 


HD 149914 


B9.58 


165 


6 


X 31.5 


16 X 61 


1 


rll000832 


3470 


HD 150193 


AO* 


150 


4 


X 6.3 


4 X 6.3 


1 


rll006208 


3470 


VV Ser 




415 


2 


X 31.5 


2 X 61 


2 


r5651200 


c2d 


LkHa 348 


Blio 


415 


4 


X 6.3 


4 X 14.7 


2 


r9831424 


c2d 


HD 163296 


A0» 


122 


4 


X 6.3 


4 X 14.7 


2 


r5650944 


c2d 


HD 179218 


B98 


244 


4 


X 6.3 


4 X 6.3 


1 


rl 1006976 


3470 


HD 190073 


AO* 


767 


4 


X 6.3 


8 X 6.3 


1 


rll007232 


3470 


LkHa 224 


A4-F9I''' 


980 


8 


X 6.3 


8 X 14.7 


2 


rl6827648 


c2d 



Note. 
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Type of background observation: 1 - dedicated background observation obtained at the same time at the on-source 
observation. 2 - "best effort" archival background observation. 
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10000 r 




-20 -10 10 20 

Ka-Kc 

Fig. 3. — Energy diagram of the rotational ladder of Hj^O. The 
numbers indicate the J rotational quantum number. Note how 
the energies of the levels generally increase within a J ladder with 
Ka — Kc- Green circles indicate the dominant upper levels of the 
13.3 and 29.5 /im line complexes. Blue stars indicate the dominant 
upper levels of the 15.17 and 17.22 /im complexes. The red square 
is the dominant upper level of the 29.85 fim line. 

minor exceptions. This is illustrated in Figure [3l which 
shows the energy levels as a function of Ka — Kc- 

Figure m shows the Spitzer SH/LH spectrum of RNO 
90, a strong water emission source, with a selection of the 
strongest line complexes identified. The energy traced by 
each transition can be easily determined by inspection of 
Figure [21 As can be seen in Figure HI most "lines" are in 
fact blends of 2-5 pure rotational water transitions, usu- 
ally with a relatively wide ran ge of excitation energies 
(see also iMeiierink et al.l I2009D . The typical structure 
of a water line complex is illustrated in Figure [5] The 
identifications include transitions that dominate the to- 
tal line blend flux, weaker lines are excluded for clarity. 
The Spitzer range includes transitions from J = 5 to 
J = 18; higher J's are present, but at a weaker level. 
Many "ortho/para" transition pairs, namely those that 
come from upper levels with Ka close to J, tend to lie 
very close in frequency, causing them to be blended at 
R=600. One example of this is the ortho lOioi Q90 
and para lOioo ^891 transitions at 20.97/im. This 
illustrates the importance of future high spectral resolu- 
tion observations in determining the ortho/para ratio(s) 
in disks. 

As much as 30% of the wavelength space is blan- 
keted in water lines above the 1% line-to-continuum level 
(jPontoppidan et 311120091) . An important consequence of 
this is that peaks tend to appear in the Spitzer spectrum 
at places where lines crowd together, and not necessarily 
where a single, particularly strong line is present. RNO 
90 is chosen as an illustrative example because it has 
one of the highest line-to-continuum ratios in the Spitzer 
sample, but it may not have particularly high contrast 
lines at higher resolving power. 

In addition to water, lines due to other molecular 
species are present in many Spitzer spectra, includ- 
ing that of RNO 90. OH exhibits a characteristic 
doublet pattern, and at least three sets of OH lines 
are visible around 23.2, 27.5 and 30.5 fj,m. The Q 
branches of C2H2, HCN and CO2 are seen at 13.7, 14.0 
and 1 4. 95 /xm, respect i vely, a s als o noted in iSalyk et al.l 
(|2008f) . lCarr fc Naiital (f200l and lPascucci et all (|2009i r 
Additional features of note are the [Ne II] line at 



12.814/zm, the H2 S(2) line at 12.279 /um, as weU as sev- 
eral H I lines across the Spitzer range, the brightest one 
being the (7-6) line at 12.371 /xm. In general, all of these 
transitions are blended with water lines, which should 
be taken into account when measuring line strengths in 
water-rich spectra. 

The water and organics line emission is spatially unre- 
solved, in contrast to the H2 lines which are sometimes 
extended (|Lahuis et al.l [20071) . This is confirmed both 
by inspection of the spectral images, as well as from the 
dedicated background observations, where available. 

4. IDENTIFICATION OF MOLECULAR TRACERS 
4.1. Detection criteria 

Due to the large number of lines, almost all of which 
are blended with other lines of similar strength, the anal- 
ysis of even a single Spitzer spectrum can be a daunt- 
ing task. Further, for sources of somewhat lower quality 
(signal-to- noise ^ 100), water emission close to the de- 
tection limit may be difficult to distinguish from noise 
and residual fringe patterns. Therefore, for the purposes 
of: 1) separating disks where water emission is clearly 
detected from non-detections and borderline cases, and 
2) defining a simple excitation temperature parameter, 
a few strong lines were selected that could serve as trac- 
ers. Such tracer line complexes were selected among fea- 
tures that are relatively isolated in the spectrum, allow- 
ing a continuum to be fitted. Experiments were made in 
which a generic LTE model for water (see paper II) was 
correlated across a wide wavelength range of the spec- 
trum. However, it was found that this did not produce 
results much different from the approach of using a few 
line complexes. The water line complexes at 15.17 and 
17.22 /im are more isolated than most, and while they 
are both blends of several lines, the excitation tempera- 
ture of the transitions contributing to each complex fall 
within a narrow range. The 15.17/xm complex traces ro- 
tational J quantum numbers of 8-10 with excitation en- 
ergies of 2300-2700 K, while the 17.22 /im complex traces 
Js of 9-12, corresponding to slightly higher energies of 
2300-3000 K. Thus, these two complexes trace roughly 
the same gas, and they are expected to correlate. 

A detection of water emission is defined as a detection 
of both 15.17 and 17.22 /im complexes at the 3.5cr level. 
In a few cases, the 15.17 /im line was not formally de- 
tected, but other water lines were strong enough to still 
warrant a clear detection. Conversely, in some cases lines 
were formally detected, suggesting the presence of water 
emission that a visual inspection could not confirm. For 
instance, a lack of other water lines in the same spec- 
tral range would lead to the source being fiagged as a 
non-detection. Due to the nature of the data set, such 
subjective analysis of a few borderline cases is unfortu- 
nately inevitable. We aimed to be conservative to en- 
sure that there are no false positives in the sample, at 
the expense of rejecting a few real detections. Hence, 
all detection rates can be considered lower limits. Cases 
that were subjectively scrutinized but where we could 
not unambigously confirm detections of, in particular, 
water are flagged in Tables [3] and [D Selected regions, 
including those containing the water tracer lines of the 
Spitzer-IRS spectra of T Tauri disks for which water is 
detected are shown in Figures [6| and [3 for the SH and LH 
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Fig. 4. — Representative identification, using the spectrum of RNO 90, of selected major molecular line complexes in the Spitzer high 
resolution spectral range (many features are left unlabeled for clarity). Unless otherwise noted, the transitions refer to the rotational 
quantum numbers Jxa Kc ™ the ground vibrational state of tlg'^O. Rotational transitions in the vibrational ground state where Ka + Kc 
is odd have ortho nuclear spin functions. 
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Fig. 5. — LTE disk model showing how the 17.22 /im H2O line 
complex breaks up into many individual lines at higher spectral 
resolution. Model spectra with resolving powers relevant for both 
Spitzer-IRS {R = 600) and JWST-MIRI (R=3000) are shown. 
Note how JWST-MIRI is expected to resolve the water line com- 
plexes (but not the individual lines). The lower curve is the spec- 
trum of TW Cha, shown for comparison. The model is taken from 
IPontoppidan et al.l l[2009h . 

modules, respectively. For comparison, spectral regions 
including the water tracer lines as well as the HCN and 
C2H2 Q branches for the Herbig Ae/Be disk sample are 
shown in Figure |S1 although these spectra are all non- 
detections. In practice, the detection rate does not de- 
pend on the exact choice of tracer line complexes, except 
in a few borderline cases. Detections of OH are based 
on the doublets around 23.0, 27.6 and 30.5 /im, span- 
ning upper level energies of ^2400-4000 K. The organics 
HCN/C2H2 and CO2 are detected via their characteristic 
Q branches. For the latter species, it is noted that there 
may be a bias against their unambiguous detection in 
sources showing strong water emission, due to blending 
and confusion. In addition, these species are identified 
based on a single blended feature making it difficult to 
rule out the occasional spurious detection, especially if 
systematic errors or residual fringing exceed calculated 
errors. A few such spurious detections were eliminated 
with visual inspection. 

The general low contrast of the molecular lines at 
R=600 relative to the strong continuum also affects the 
detection rate. Specifically, the fringe and flat field resid- 
uals, in addition to the photon statistics, scale with the 
continuum, so even for very high theoretical signal-to- 
noise ratios, there is a limit to how low the line-to- 



continuum ratio can be to allow the detection of lines. 
For our data, this "fidelity" limit, i.e. the line-to- 
continuum contrast where the data systematics become 
larger than the pure photon statistics, while difficult to 
quantify, seems to be better than 0.5%, with several 3.5cr 
detections at the 1-2% level. 

Figure ^ shows the relation between the continuum 
level and the flux level of the water tracer. The ability to 
detect a line of constant flux diminishes with increasing 
flux level, indicating the "fidelity" limit. On the other 
hand, a significant number of line detections are made 
in T Tauri disks with apparent brightnesses higher than 
those of many Herbig Ae/Be disks. 

Using these criteria, water emission is detected in 22 
disks, HCN in 25 disks, C2H2 in 17 disks, CO2 in 20 
disks and OH in 18 disks, as summarized in Figure [TU] 
and Table [31 Figure [TUj also shows the detection rates 
of CO in the rovibrational fundamental band at 4.7 /im. 
The CO observations are described in Paper II. 

4.2. H2O excitation tracers 

The 15.17 and 17.22 /.tm water line complexes chosen 
to act as signposts for molecular emission trace temper- 
atures in the middle of the range of the 10-36 /xm water 
lines. However, the spectra also include line complexes 
from strong transitions with much lower excitation tem- 
peratures. Ratios of line fluxes between the high and 
low excitation line complexes are an important diagnos- 
tic of the nature of the molecular emission; indeed one 
of the first questions that we ask, based on the present 
sample, is whether the excitation conditions are simi- 
lar for all disks, or whe ther there are si g nifica nt differ- 
ences. As discussed in iMeiierink et al.l ()2009() for the 
particular case of mid-infrared water lines, transitions 
with lower excitation temperature trace larger radii of 
the disk. Hence, line ratios can be affected by both exci- 
tation conditions and the radial abundance distribution 
of the species in question. In any case, it will be nec- 
essary to identify lines with the widest possible baseline 
in excitation energies. A good low energy tracer is the 
ortho 72 5 ~^ 61 6 line at 29.85 /im with an excitation 
energy of HOOK (see Figure[3|). Spectra with higher re- 
solving power will be able to separate lines tracing levels 
at least down to J = 5; Ka — Kc = 3 and E ~ 900 K, 
but these are too blended with higher excitation lines 
at the Spitzer resolution to be clean tracers of specific 
conditions. Very high excitation tracers include lines 
at 13.3 /im (J=15-16, E = 4400 - 5200 K) and 29.5 /im 
(J=16-18, E = 4600 - 4800 K). 

Table [31 summarizes the integrated fluxes of the 15.17, 
17.22 and 29.85 fim line complexes. The line fluxes were 
calculated by fltting a Gaussian superposed on a lin- 
ear continuum to the data. The selected line complexes 
appear to be spectrally unresolved, so the centers and 
widths of the Gaussian were kept constant for each line 
complex, corresponding to the theoretical center of the 
line complex and the nominal spectral resolving power of 
Spitzer-IRS. The wavelength range of the flt (determin- 
ing where the continuum is constrained) is indicated by 
the models shown in Figures [51 and [7] 

Figure [TT] shows the distribution of detected water 
sources in a 3- line plot, matching the high excitation 
15.17 and 17.22 /im complexes with the low excitation 
29.85 /zm complex. Since the low excitation line probes 
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Fig. 6. — Selected regions of the Spitzer-IRS SH module of T Tauri stars with detected H2O. From left to right, the spectra are centered 
on: 1) the acetylene (C2H2) and HCN Q-branches, 2) the 15.17/xm H2O line complex, and 3) the 17.22 /jm H2O line complex. The lines 
are marked with a Gaussian fit (red curves). The spectra are in order of decreasing line-to-continuum ratio of the 17.22 /im complex. 
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Fig. 7. — As for Figure|6l but for Spitzer-LH data, including the 21.76, 29.85 and 33.0 /^m H2O complexes, as well as a region centered 
on the OH doublets around 30.5 fim. 
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Fig. 8. — Selected Spitzer-IRS spectral regions from the full HAe/Be disk sample for comparison with Figures [6] and [7] The spectra are 
ordered according to signal-to-noise, with the lowest signal-to-noise at the top. The dashed lines indicate the locations of the tracer water 
line complexes, except where otherwise noted. 



the disk surface at larger radii than the high excitation 
hnes, the location of a source in this diagram could pro- 
vide an indication of the radial water abundance struc- 
ture of the disk surface: the smaller the 29.85/17.22 ra- 
tio, the sharper the cutoff of water em ission beyond a 
certai n radius. This was discussed in iMeiierink et al.l 
(|2009() . who presented generic non-LTE models of emis- 
sion spectra, one for a constant abundance of water, and 
one in which the surface water abundance was lowered by 
6 orders of magnitude beyond a radius of 0.7 AU, corre- 
sponding to the location of the mid-plane snow-line. The 
arrow in Figure [11] indicates the approximate range and 
direction of increasing "cold finger depletion" of the sur- 
face water abundance beyond the mid-plane snowline. It 
is seen that essentially all disks are located between these 
two extreme cases, with only a few disks exhibiting line 
ratios that are consistent wit h a high abundan ce of water 
throughout the disk surface. IMeiierink et al.l |2009) sug- 
gests a mod e l, bas ed on a few H2O Spi t zer sp ectra from 
ISalvk et al.l (12001 and ICarr fc Naiital (l2008l) . in which 



the disk surface water vapor may be strongly depleted 
at radii larger than that of the midplane snow line due 
to vertical transport of water. Purely chemical effects 
may also contribute to generate structure in the radial 
surface abundance structure. In any case, the data, as 
shown in Figure [TTl indicate that line ratios cluster for 
many sources, but significant variation does exist in the 
sample, with one possible interpretation being that this 
is due to structure of the radial emission profile of water. 

5. DISCUSSION 

5.1. Cooling balance 

The temperature structure of protoplanetary disks is 
determined by a detailed balance of heating and cool- 
ing processes. In the surface layers, collisional exchange 
between the gas and the dust is not sufficient to main- 
tain a temperature equilibrium between the dust and the 
gas, leading to elevated gas temperatures. The degree to 
which the gas is superheated depends sensitively on the 
efficiency of line cooling. As a consequence, the finding 
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TABLE 3 

Line fluxes and molecular detections from the T Tauri star sample 



Source name 


15.17 ;im=' 


17.22 lira 


29.85 lira 




H2 0'^ 


OH 


HCN 


C2H2 


C02 




lO""*^^ erg cm~^ s""*^ 


10""*^^ erg cm~^ s""*^ 


10""^^ erg cm~^ 


10-^Lq 












LkHa 270 


< 0.36 


< 0.34 


0.82 ±0.07 


< 1.2 


(0) 











1 


LkHa 271 


< 0.15 


< 0.33 


< 0.68 


< 1.0 

















LkHa 326 


0.47 ±0.05 


0.62 ±0.05 


0.30 ±0.04 


1.3 


1 


1 


1 





1 


LkHa 327 


1.00 ±0.27 


2.43 ±0.25 


2.31 ±0.15 


5.9 


1 


1 








1 


LkHa 330 


< 0.23 


< 0.29 


< 0.98 


< 1.3 

















LkCa 8 


< 0.17 


0.27 ±0.05 


0.71 ±0.06 


< 0.3 


(0) 














IQ Tau 


0.57 ±0.10 


0.90 ±0.09 


< 0.51 


0.8 


1 





1 


1 





V710 Tau 


< 0.17 


< 0.32 


< 0.69 


< 0.4 








1 


1 


1 


AA Tau 


0.48 ±0.06 


1.48 ±0.06 


0.74 ±0.03 


0.9 


1 


1 


1 


1 


1 


CoKu Tau/4 


< 0.29 


< 0.33 


< 0.39 


< 0.3 

















DN Tau 


< 0.27 


< 0.27 


< 0.59 


< 0.3 

















FX Tau 


< 0.34 


< 0.43 


< 0.75 


< 0.4 








1 








DR Tau 


4.53 ±0.19 


7.13 ±0.19 


3.73 ±0.10 


4.4 


1 


1 


1 


1 


1 


SX Cha 


0.78 ± 0.06 


1.21 ±0.06 


1.01 ±0.04 


1.3 


1 


1 


1 








SY Cha 


< 0.09 


0.12 ±0.03 


0.43 ±0.02 


< 0.2 


(0) 





1 


1 


1 


TW Cha 


0.61 ±0.02 


1.04 ±0.03 


0.55 ±0.02 


1.2 


1 


1 


1 


1 





VW Cha 


1.99 ±0.07 


3.54 ±0.08 


1.81 ±0.04 


3.6 


1 


1 


1 


1 


1 


VZ Cha 


0.94 ±0.05 


1.39 ±0.05 


0.75 ±0.02 


1.4 


1 


1 


1 


1 





WX Cha 


1.09 ±0.04 


1.61 ±0.04 


0.90 ±0.02 


1.6 


1 


1 


1 


1 


1 


XX Cha 


0.40 ± 0.03 


0.46 ± 0.03 


0.40 ±0.02 


0.5 


1 


1 


1 


1 


1 


T Cha 


< 0.32 


< 0.31 


< 0.92 


< 0.08 

















Sz 50 


< 0.12 


< 0.13 


< 0.18 


< 0.2 














1 


HD 135344B 


< 0.32 


< 0.35 


< 2.11 


< 0.2 

















HT Lup 


< 1.75 


< 1.72 


1.17 ±0.15 


< 1.1 


(0) 











1 


GW Lup 


< 0.17 


< 0.18 


< 0.34 


< 0.3 














1 


GQ Lup 


0.71 ±0.06 


1.39 ±0.06 


1.33 ±0.04 


1.3 


1 


1 


1 








IM Lup 


< 0.24 


< 0.24 


< 0.30 


< 0.2 














1 


HD 142527 


< 1.96 


< 2.22 


< 5.79 


< 0.5 

















RU Lup 


2.50 ±0.13 


3.49 ±0.14 


2.41 ±0.08 


2.8 


1 





1 





1 


RY Lup 


< 0.65 


< 0.74 


< 0.97 


< 0.7 

















EX Lup 


< 0.37 


1.28 ±0.12 


4.05 ±0.07 


2.4 


1 


1 











AS 205 


11.55 ±0.65 


18.19 ±0.74 


9.35 ±0.48 


8.9 


1 


1 


1 


1 


1 


Haro 1-1 


< 0.23 


< 0.22 


< 0.67 


< 0.3 

















Haro 1-4 


< 0.38 


< 0.47 


< 0.75 


< 0.4 








1 


1 





VSSGl 


1.54 ±0.18 


2.63 ±0.17 


2.78 ±0.19 


1.9 


1 





1 


1 


1 


DoAr 24E 


2.65 ±0.11 


4.42 ±0.12 


2.54 ±0.13 


2.5 


1 


1 


1 


1 


1 


DoAr 25 


< 0.21 


< 0.51 


< 0.53 


< 0.3 








1 


1 





SR 21 


< 1.04 


< 1.29 


< 3.90 


< 0.2 

















SR 9 


< 1.05 


< 0.70 


1.32 ±0.08 


< 0.5 


(0) 














V853 Oph 


< 0.28 


0.55 ±0.08 


1.14 ±0.05 


0.6 


1 


1 


1 


1 





ROX 42C 


< 0.38 


< 0.37 


< 0.68 


< 0.3 

















ROX 43A 


< 0.85 


< 0.92 


< 1.36 


< 0.7 

















Haro 1-16 


0.83 ±0.07 


1.56 ±0.08 


1.44 ±0.06 


1.0 


1 


1 


1 








Haro 1-17 


< 0.14 


< 0.14 


< 0.35 


< 0.2 

















RNO 90 


5.83 ±0.24 


10.10 ±0.25 


5.86 ±0.14 


5.8 


1 


1 


1 


1 


1 


Wa Oph/6 


1.57 ±0.07 


1.54 ±0.07 


1.06 ±0.05 


0.8 


1 


1 


1 





1 


V11210ph 


1.12 ±0.30 


2.69 ±0.32 


2.54 ±0.17 


2.3 


1 














EC 82 


< 0.31 


< 0.35 


< 0.65 


< 4.9 


















^ All upper limits are 3.5 cr; errors arc la. 

^ The integrated line flux due to H2O in the Spitzer wavelength range. Water has a multitude of lines outside the observed speetral range, 
ineluding the strong rovibrational band around 6 ^m. Extrapolating to all water lines may signifleantly increase the total cooling rate. 

(0) indicates a tentative detection of water based on an inspection by eye. Some times, but not always, these borderline detections are 
accompanied by formal detections of speciflc line complexes. Conversely, in a few cases, formal detections were made that could not be 
confirmed by a visual inspection, likely due to interference from data artifacts. 



that the mid-infrared wavelength range in typical proto- 
planetary disks is blanketed in molecular lines will have 
important consequences for models of disk structures. 

Current disk models find water cooling rates of a few 
x10~^Lq and individual cooling line luminosities from 
atomic species, as wel l as rotational lines of CO and H2, 
are 10"* - IQ-'^Lq (|Gorti fc Hollenbachl[200l ). How- 
ever, we measure integrated water line luminosities in 
the 10-36 /im region that are 1-2 orders of magnitudes 
higher. Table [3] presents the 10-37 /im integrated H2O 
line luminosities for the disks in which H2O has been de- 
tected. The line lumi nosities were determine d by scaling 
the fiducial model of iMeijerink et alJ (|2009| ) to the SH 



and LH spectra, where different scaling values were al- 
lowed for the two spectral ranges. Typical values result- 
ing from this range from < 5 x 10~'' to almost 10"^ L©. 
Extrapolating to include lines outside the Spitzer range 
may increase these numbers by a factor 2 or more. Con- 
sequently, the detection of the mid-infrared molecular 
spectra from protoplanetary disks increases the total 
disk-averaged cooling rates by 1-2 orders of magnitude, 
relative to recent models, for gas temperatures below 
< 2000 K; the infrared cooling budget is dominated by 
water in T Tauri stars. It should be stressed, however, 
that most of the water lines in the Spitzer range are likely 
formed in the innermost (< a few AU) regions of the disk. 
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TABLE 4 

Line fluxes and molecular detections from the Herbig Ae/Be star sample 



Source name 


15.17 ^m'' 


17.22 fim 


29.85 fim 




H2O'' 


OH 


HON 


C2H2 


c 




















HD 36112 


< 1.58 


< 1.93 


< 3.84 


< 5.8 


(0) 














HD 244604 


< 0.86 


< 0.83 


< 0.91 


< 5.9 

















HD 36917 


< 0.94 


< 1.01 


< 1.44 


< 9.7 

















HD 37258 


< 0.80 


< 0.89 


< 0.70 


< 6.2 


(0) 














BF Ori 


< 0.62 


< 0.59 


< 0.72 


< 3.9 

















HD 37357 


< 0.78 


< 0.79 


< 1.11 


< 11.3 

















HD 37411 


< 0.38 


< 0.38 


< 1.05 


< 8.0 

















RR Tau 


< 0.69 


< 0.68 


< 0.93 


< 172 


(0) 














HD 37806 


< 2.01 


< 2.09 


< 2.45 


< 38 


(0) 














HD 38087 


< 0.22 


< 0.27 


< 0.87 


< 5.7 

















HD 38120 


< 2.15 


< 2.62 


< 4.39 


< 58 

















HD 50138 


< 4.57 


< 4.99 


7.81 ±0.89 


< 84 


(0) 














HD 72106 


< 0.86 


< 0.88 


< 1.16 


< 4.0 


(0) 














HD 95881 


< 1.78 


< 1.75 


< 2.37 


< 1.4 

















HD 98922 


< 4.66 


< 4.83 


< 3.96 


< 346 


(0) 














HD 101412 


< 0.86 


< 0.82 


< 0.82 


< 1.1 


(0) 











1 


HD 144668 


< 2.42 


< 2.45 


4.25 ± 0.53 


< 8.7 


(0) 














HD 149914 


< 0.21 


< 0.16 


< 0.36 


< 0.16 

















HD 150193 


< 2.35 


< 2.61 


< 4.04 


< 5.8 


(0) 














VV Ser 


< 0.94 


< 1.15 


0.99 ±0.14 


< 11.7 


(0) 


(0) 











LkHa 348 


< 1.89 


< 1.78 


< 1.98 


< 13.3 

















HD 163293 


< 2.42 


< 2.73 


4.34 ± 0.39 


< 4.3 


(0) 














HD 179218 


< 3.14 


< 3.76 


< 6.26 


< 54 

















HD 190073 


< 1.64 


< 1.66 


< 1.75 


< 65 

















LkHa 224 


< 1.52 


< 1.57 


< 2.04 


< 219 


















^ All upper limits arc 3.5 cr; errors arc 1 a. 

^ (0) indicates a tentative detection of water based on an inspection by eye. Some times, but not always, these borderline detections are 
accompanied by formal detections of specific line complexes. Conversely, in a few cases, formal detections were made that could not be 
confirmed by a visual inspection, likely due to interference from data artifacts. 
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Fig. 9. — Relationship between the 17.22 fim integrated line flux 
and the continuum level. The filled symbols show detections, while 
empty symbols with arrows indicate upper limits. Squares, stars 
and circles represent T Tauri stars, transitional objects and Herbig 
Ae/Be stars, respectively. 

the outer regions of which may be less affected by wa- 
ter as an exceptionally strong coolant. An immediate 
consequence of more efficient cooling seems to be that 
the surface layers of the disks will have a larger column 
of relatively cooler gas, facilitating the survival, and ob- 
servability, of a rich chemistry in the disk surface. 




F G 
Spectral Type 

Fig. 10. — Detection rate as a function of spectral type for the 
strongest infrared molecular tracers, with the number of each spec- 
tral type shown above. The detection rate drops dramatically for 
spectral types earlier than G/F, except for CO. 

5.2. Lack of detectable molecular emission from Herbig 
Ae/Be stars 

Arguably the strongest observation, based on the cur- 
rent sample, is that the detection rate of molecular emis- 
sion from Herbig Ae/Be stars with Spitzer-IRS is appar- 
ently very low. The upper limits for detections of the 
water tracers are given in Table |31 In fact, out of 25 
Herbig Ae stars, there is not a single detection of water, 
HCN, C2H2 or OH, given our detection criteria; CO2, 
which is detected in HD101412, provides the sole excep- 
tion. In comparison, the detection rate for disks around 
stars with spectral types later than F is '~40%. Could 
this effect be due to some bias? First of all, as Figure [1] 
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V 



15.17/17.22 

Fig. 11. — 3-line plot of the observed sample using two high exci- 
tation line complexes (15.17 and 17.22 /xm) and one low excitation 
transition (29.85 /^m). The high excitation lines probe the inner- 
most radii of the disk with temperatures of T = 1000 — 2000 K, 
while the lower excitation line probes larger radii with T = 
300 — 500 K. The lower part of the plot signifies low water abun- 
dances at larger radii, while the upper part of the plot indicates a 
constant surface abundance. 




T Tauri star WX Cha X 1.0 




Wavelength [micron] 

Fig. 12. — Comparison of typical spectra from three classes of 
objects: Transitional disks, Herbig Ae/Be stars and classical T 
Tauri stars. The spectra have been scaled for clarity (not shifted), 
preserving the line-to-continuum ratios. 

shows, the sample of Herbig stars are somewhat brighter 
than the T Tauri stars, by up to an order of magnitude. 
This indicates that the signal-to-noise ratios of the Her- 
big star spectra are unlikely to be significantly lower than 
those of the T Tauri stars. Figure fOl shows a comparison 
between the full IRS spectra of a Herbig star, a transi- 
tional disk and two classical T Tauri stars. Further, in 
Figure [13] a comparison of the line-to-continuum ratios 
of the H2O line tracers at 15.17 and 17.22 /im is given 
for the full sample. It is seen than the Herbig stars have 
3.5(7 upper limits on their line-to-continuum ratios that 
are systematically smaller by a factor of 5-10 than the 
ratios in T Tauri stars where water is detected. 

It is interesting to note that some Herbig Ae/Be disks 
do show tentative low-level emission features that may be 



. Sp. type =< F (CTTS) 

Sp. type =< F (CTTS) upper limits 
^ Sp. type >- A (HAeBe) upper limits 



— ^2r' 



£. 0.01 : 



0.01 0.10 

peak 15.17/continuum 

Fig. 13. — Distribution of line-to-continuum ratios of the 15.17 
and 17.22 fim line complexes. The sample has been split into T 
Tauri stars (spectral types later or equal to F) and Herbig Ae stars 
(spectral types earlier than F). 

due to water and OH at longer wavelengths, specifically 
in the LH module. Such tentative detections made using 
the 24.9-25.5 /im OH/H2O complex are displayed in Fig- 
ure [14] and compared to the water spectrum from TW 
Cha. We describe them as tentative because they do not 
unambiguously match a water model over a wider range 
of wavelengths, although that could be explained by noise 
or systematics that vary with wavelength or spectral or- 
der. The 25 /im complex covers the X7r3/2 ^'^3/2 10-5 
and X7r2/2 ~^ '^''^1/2 9-5 OH lines and water lines span- 
ning excitation energies of 1500-2500 K, or somewhat 
cooler than those traced by the 15.17 and 17.22 /im com- 
plexes. If real, these detections still represent line-to- 
continuum ratios significantly lower than those of T Tauri 
disks, but suggest that, with a modest improvement in 
data quality, molecular emission features should also be 
detected in Herbig Ae/Be disks. Given that the tenta- 
tive detections are at the same level of the current data 
systematics, it is not possible to further analyze them. 
Their presence, however, do suggest that the optically 
thick H2O lines at far-infrared wavelengths may be de- 
tected in Herbig Ae/Be disks by Herschel. 

While the mid-infrared molecular line-to-continuum 
ratios are much smaller in Herbig Ae/Be stars relative 
to those in T Tauri stars, is it possible that the line 
fluxes are similar, but veiled by the stronger continuum 
fluxes of the Herbig stars? Most disks around Herbig 
Ae/Be stars, as well most transitional disks, do show 
strong lines from CO at the rovibrational band at 4.7 /im 
(jBlake fc Boogert|[200l iSalvk et al.|[200l and Paper II). 
New ground-based CO observations of the full sample 
are discussed in greater detail in Paper II. 

A comparison of line- to-continuum ratios for 
CO M-band li n es fro m INaiita et all ([2003) and 
iBlake fc BoogertI (|2004[) demonstrates lower line-to- 
continuum contrast, on average, for Herbig Ae/Be vs. 
T Tauri disks. The reason for this difference is still 
unclear, but if the water and other molecules behave 
similarly, this may explain the general lack of high 
contrast molecular line emission from Herbig Ae/Be 
stars. However, the difference for CO is only a factor 
of 2 — significantly smaller than the minimum difference 
required for water. 

It is therefore tempting to conclude that the observed 
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24.0 24.5 25.0 25.5 26.0 
Wavelength [micron] 

Fig. 14. — Tentative detections of line complexes due to a combi- 
nation of H2O and OH lines in the sample of Herbig Ae/Be disks, 
specifically, the complex located between 24.9 and 25.5 fira. The 
top spectrum of TW Cha is displayed for comparison with a source 
with clear detections. The two vertical lines indicate the locations 
of two OH doublets. Due to the low line-to-continuum of the lines 
and the possibility of interference from systematics at the same 
level, we refrain from analyzing these further. 

lower line-to-continuum ratios in Herbig Ae/Be disks are 
produced by lower molecular abundances, although other 
alternative explanations are possible. Differences in dust 
properties could affect line strengths, or the heating rate 
of the molecular layer may not scale linearly with stellar 
luminosity. Here, we discuss these possibilities. 

5.2.1. Photochemtstry 

In principle, the harsher radiation field from A and B 
stars effectively photodissociates molecules in unshielded 
regions of the disk surface at 1-2 AU, which could result 



in lower line fluxes relative to the disk luminosity. If pho- 
tochemistry is important, one would expect to see strong 
chemical changes with spectral type, such as an increas- 
ing OH/H2O abundance ratio for disks around earlier 
type stars. Some early indications of this have already 
been found in two studies reporting detections of O H, but 
not H 2O, for one Herbig Ae/Be star by Mandell et al.l 
(|2008[) and in a transitional disk (jNaiita et al.„2010t) . A 
dependence of inner disk emission on spectral type be- 
tween late type T Tauri stars and brown dwarfs has also 
been found for the HCN/C2H2 ratio by Pascucci et al] 
(j2009( ). although in this case photochemistry is unlikely 
to be the cause of the observed difference. Photochem- 
istry may also give rise to a dichotomy between molecules 
(like CO) that can efficiently self-shield and those that 
cannot; molecular self shielding may be a very strong ef- 
fect in disks, i.e. even stronger than in molecular clouds 
because the continuum dust opacity is expected to be 
small due to a depletion of dust opacity in the disk sur- 
face following dust growth and settling to the midplane. 

Water is an interesting case for self-shielding that was 
recently studied for the case of protoplanetary disks 
by BethcU & Bcrgin (2009). These authors show that 
H2O, as well as OH, can in fact self-shield, although 
at higher column densities than CO. Specifically, wa- 
ter self-shields at column densities of A^HaO 2 x 
10^^ cm-^ while CO se l f-shiel d s at Ngn ~ 5 x 10^ ^ cm-^ 
()van Dishoeck fc Black! 119881 : iVisser et al.] I2009D in a 
dust-free environment, and assuming abundances of ~ 
10~* relative to H. This suggests that CO may survive 
closer to the star in regions unshielded by dust, pre- 
dicting that low excitation water lines may still be seen 
from larger radii in Herbig disks, even if water is absent 
from warm gas. CO emission is generally observed in 
the inner regions of cl assical Herbig Ae/Be disks (e.g. 
iBlake fc Boogerj l2004f ). including the majority of the 
disks in the sample presented here (Paper II) . Note that 
some disks around Herbig Ae/Be stars are also part of 
the class of transitional disks, and some of those appear 
to ha ve strongly depleted CO abundances in their inner 
disks ()van der Plas et al.ll2009HBrittain et al.ll2009D . but 
that may be related to their evolutionary stage, rather 
t han to an intr i nsic pr operty of Herbig Ae/Be disks. 

iWoitke et al.l ()2009f ) modeled the abundance of wa- 
ter vapor throughout the disk surrounding a Herbig Ae 
star, including a detailed treatment of UV-driven pho- 
tochemistry, and found that the inner disk abundance 
is high, if shielded, reaching 10^^ — 10^"* per hydro- 
gen nucleus, and can maintain abundances of 10^^H~^ 
in higher, unshielded layers. While there are no pre- 
dicted line fiuxes for the Spitzer wavelength range, this 
model does produce strong lines at wavelengths > 70 fim. 
Line spectra in the Spitz er w indows were generated by 
iPontoppidan et al.l (|2009( ) and lMeiierink et al.l (|2009D for 
similar water abundances, but for a lower mass star, 
which resulted in strong water lines. A detailed radiative 
transfer comparison of mid-infrared molecular lines from 
the disks around stars of varying stellar mass is therefore 
needed, and current chemical disk models may have to 
consider additional effects to explain the observed differ- 
ence between T Tauri and Herbig Ae/Be disks. 

5.2.2. Relative scaling of line and continuum luminosities 
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As is seen in Tables |3] and 21 the upper limits on the to- 
tal H2O line luminosities for the Herbig star sample tend 
to be higher than the line luminosities for the T Tauri 
stars with detected lines, with a few exceptions. Because 
the mid-IR continuum level is generally proportional to 
the stellar luminosity, it is possible for the molecular 
spectra to remain undetected in the Herbig Ae/Be stars 
if the molecular line luminosity is a significantly weaker 
function of stellar luminosity. Ultimately, the scaling of 
line fluxes with the stellar luminosity may be related to 
how the regions of the disk surface that form the molecu- 
lar lines are heated. Three principal sources of disk sur- 
face heating have been identified: Heating by the stel- 
lar optical/infrared continuum of dust and subsequent 
coUisional coupling to the gas, photo-electric heat ing by 
UV p hotons Q onkhcid ct al. 2004; Kamp & DuU emondl 



l2004t iNpm ura &: Millar 2005) and ionization heating by 



X-rays (|Glassgold et al. 2004). 

The dust continuum is generated by the first of these 
processes, and if that dominates the heating of the line- 
forming layer, the line luminosity can be expected to 
scale similarly with stellar luminosity, resulting in no 
significant line-to-continuum dependence with spectral 
type. The strength of the UV radiation field is a strongly 
increasing function of stellar effective temperature (and 
for T Tauri stars, of accretion rate) and would, in the 
absence of chemical differences, be expected to increase 
line-to-continuum ratios with increasing stellar effective 
temperature. The X-ray luminosity, on the other hand, 
is a much weaker function of stellar type, although highly 
variable, and may, in median, only differ by 1-2 orders 
of ma gnitude between T T auri stars and Herbig Ae/Be 
stars (jFeigelson et al.l [20051 ). If the molecular line lumi- 
nosity is required to be, at most, constant with spectral 
type, and in the absence of differences in chemistry and 
excitation. X-ray heating of the molecular layer would 
likely be the best candidate for determining the line lumi- 
nosity. Regardless, photo-electric heating becomes weak 
at column densities in excess of 10^^ cm~^ and tends 
to be a ssociated with layers w ith low molecular abun- 
dances (jDuLlemond et ani2007l ). Given that the column 
densities of the layers emitting the mid-i nfrared molec- 
ular l ines are likely hig her > 10^^ cm~^ (jCarr fc Naiital 
I2OO8I: ISalvk et aLll2008[ ). the mechanism for gas heating 
in the uppermost layers may not be important. If photo- 
electric heating, by X-rays or otherwise, is not the domi- 
nant heating mechanism in the line-forming regions, dust 
coupling is the main candidate, leading to a strong de- 
pendence of the line luminosity on stellar spectral type. 
In conclusion, we find it unlikely that the molecular line 
luminosity is a sufficiently weak function of stellar lu- 
minosity, in the absence of chemical or other systematic 
differences in the structures of disks around T Tauri rel- 
ative to those surrounding Herbig Ae/Be stars. 

5.2.3. Dust properties 

The lines from Herbig stars could be intrinsically 
weaker due to a relative difference in the gas-to-small- 
dust ratio in the disk atmosphere leading to a reduced 
column of observable molecules. This would require dif- 
ferent relative rates of vertical mixing and settling of 
grains for the two types of disks. Specifically, if the disk 
surfaces of Herbig Ae/Be stars are more turbulent, the 
water column above the r = 1 surface of dust may not be 



sufficien t to form strong lines. This scenario was mod- 
eled bv iMeiierink et aiT()2009D . who indeed found that 
the observed line strength in T Tauri disks required high 
gas-to-dust ratios of ^ 10"*. If the gas and dust in Herbig 
Ae/Be disks are mixed closer to the canonical gas-to- 
dust ratio of 100, this would translate to suppressed line 
strengths. The higher ionizing radiation fields around 
Herbig Ae/Be stars may allow for a better coupling of 
the disk surface to the stellar magnetic field, thus ac- 
tivating the magneto-rotational-instabil ity and drivin; 
a larg er degree of surface turbulence (B albus fc Hawle; 
Il998| ) . A highly turbulent disk would also work to coun- 
teract the effect in which water (and therefore oxygen) 
is deplet ed from the disk surf ace due to freeze-out and 
setthng C Meiierink et al.ll2009D . 

A difficulty with this scenario is that it predicts that 
lower line-to-continuum ratios should be observed across 
the infrared range for all observed molecular lines, includ- 
ing the CO ro- vibrational lines at 4.7 /xm. Since these 
lines are generally observed in Herbig Ae/Be disks, with 
colum n densities similar to those observed in T Tauri 
disks (iNaiita et al.l2003HBlake fc Boogertl2004D . it is not 
obvious that there is a general difference in gas-to-dust 
ratio in T Tauri disks relative to Herbig Ae/Be disks. 
Also, X-rays, as compared to UV, are more effective in 
generating ionization in la yers deep enough to fac ilitate 
efficient turbulent mixing (|Igea fc Glassgoldl 119990 . but 
X-rays will not scale as rapidly with stellar type as UV 
as discussed section 15.2.21 arguing against this mecha- 
nism. 

5.3. Lack of molecular emission from transitional disks 

Transitional disks are protoplanetary disks with a 
deficit of infrared emission due to a paucity of dust in 
the inner regions (Strom et al.l ll989l). This class of disks 
is thought to be in the process of clearing out their in- 
ner regions through various processes, possibly including 
planet formation. There are 5 transitional disks in this 
sample, LkHa 330, SR 21, T Cha, HD135344B and CoKu 
Tau/4. However, none show H2O, HCN, C2H2 or CO2 
emission at Spitzer wavelengths detectable in our data. 
In contrast, LkHa 330, SR 21 and HD135344B show 
CO in emi ssion in the fundamental rovibrational ban d 
at 4.7 /xm (jPontoppidan et al.l [20Q1 ISalvk et al.l [20091) . 
Given the detection rates of molecular emission around 
regular disks, the current sample of transitional disks 
is too small to draw any firm statistical conclusions. 
HD135344B is an F star, so may be of sufficiently early 
type to lack molecular emission for the same reasons that 
the A and B stars lack it. CoKu Tau/4 is not a true disk 
in transition since dust and gas appears to have been 
cleared out by a s tellar companion and not as a result 
of disk evolution (jlreland fc Kraud 120081 ). The chance 
that the remaining 3 transitional disks do not show H2O 
emission by coincidence, given the detection rate of the 
remaining disks, is only a few percent, but we do not con- 
sider this enough to conclude that no transitional disks 
will show molecular emission (especially if deeper spectra 
are acquired). 

The absence of strong molecular emission from these 
systems is similar to the results of high signal-to-noise 
IRS spectroscopy of the transitional object TW Hya 
(|Naiita et al.ll2010| ). The spectrum shows a striking lack 
of strong emission features of H2O, C2H2, and HCN in 
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the 10-20 fim region, although weak emission is detected 
from other molecules (H2, OH, CO2, HCO"'", and ten- 
tatively CH3). Najita et al. describe how the lack of 
strong molecular emission is consistent with the possibil- 
ity that the inner disk has been cleared by an orbiting 
giant planet, although chemical and/or excitation effects 
may be responsible instead. 

A lack of strong molecular emission from transitional 
disks is consistent with significant photodestruction due 
to insufficient shielding from dust in the optically thin 
inner region of the disks. In this scenario, CO is still 
seen in rovibrational transitions due to a combination of 
survival through self-shielding and the fact that CO has a 
mechanism for UV fluorescent excitation that is efficient 
also for cold, low density gas. The latter could be the 
case for SR 21, for which the CO gas originates from a 
ring at 6-7 AU, distances at which UV fluorescence of CO 
is likely operating. The CO gas from HD 135344B, on the 
other hand, originates m ostly from small radii (< 1 AU 
iPontoppidan et al.|[2008[) . 

5.4. CO2 sources 

A significant sub-class of disks has been identified, con- 
sisting of sources showing strong 14.98 fim CO2 emission 
from the Q-branch of the fundamental bending mode, 
V — (0,1,0) — > (0,0,0), but which show no other de- 
tectable molecular emission in the Spitzer wavelength 
range. Because the identification with CO2 is based on 
only one feature, it was confirmed that the emission is 
present in both the A and B nod positions and that 
it is not due to an artifact from the edge of order 13, 
which begins at 15.08 ^m. Some oxygen-rich asymp- 
totic gia nt branch (AGB) stars also show strong CO2 
emission (lJusttanont et al.lll998HCami et al.l l2000). but 
with much higher optical depth leading to strong features 
from combination bands at 13.87 and 16.18/zm. These 
(intrinsically weaker) bands are absent from protoplane- 
tary disk spectra, suggesting relatively low optical depth 
of the bending mode. The data from the CO2 disks are 
shown in Fig. I15[ and compared to a generic disk model 
from (|Meiierink et al.i i2003) generated using the radia- 
tive transfer code RADLite (jPontoppidan et all I2009D . 
assuming level populations in local thermodynamic equi- 
librium and a CO2 abundance of 6 x 10~^ relative to H2. 
This abundance is 1-2 orders of mag nitude lower than 
that inferred from chemical models (jWillacv fc Woodi 
[2009I) . but given that the critical density for the CO2 
bending mode up per state is high (10^° — lO^^cm"^, 
ICastle et all l2006h non-LTE calculations will probably 
result in higher abundances. Spectra at higher resolution 
should reveal strong R and P branches in the 14-16 /^m 
range. Further chemical and radiative transfer modeling 
is required, not so much to explain the presence of the 
CO2 emission, but why that from water, HCN and C2H2 
is absent. Are the abundances of these molecules really 
lower relative to CO2, or is the absence of emission due 
to small differences in the density and temperature struc- 
ture of the molecular surface layer of these disks causing 
other molecules to be subthermally excited or shielded 
by dust? 

6. CONCLUSIONS 

We have found that a large fraction of protoplan- 
etary disks around low-mass to Solar-type stars have 
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Fig. 15. — Detected CO2 emission lines at 14.98 /^m for sources 
with no other unambiguous molecular emission. For clarity, only 
every third error bar is plotted. A gener ic model of CO2 emission 
from a disk, generated using RADLite JPontoppidan et al.l 120091 ) 
is shown in the lower, left panel. 

mid- to far-IR spectra that are blanketed in emission 
lines from a wide ra nge of molecules at temperatures of 
500-1000 K (Fig. El ICarr fc NaiitU (jlfflSh . ISaMc et al.l 
(2008)). The lines have been shown to be excited in the 
0.1-10 AU regio n of the disk, corresponding to the vlanet- 
forminq zone (ICarr fc Naii tal [20081: ISalvk et al.l [20081: 
IPontoppidan et al.l l2009: Me iierink et al.ll2009D . The key 
conclusion is that such emission appears to be common 
in the disks around Sun-like stars. It is not a rare or ex- 
otic phenomenon, but is rather an unequivocal statement 
that the chemical environment in the planet-forming 
zone is extremely rich. While just a handful of abundant 
molecules have been identified so far due to the relatively 
low spectral resolution of space-based instruments, the 
Spitzer results demonstrate the concurrent presence of 
0-, C- and N-chemistry. With these ingredients and at 
the densities and temperatures of the inner disks, the 
chemistry must be highly complex. Given these data, 
it is the expectation that IR studies at higher spectral 
resolution will reveal many more molecular species. 

Why have these complex infrar ed emission spec- 
tra not previously been seen (e.g., IMeeus et al.l 120011 : 
IKessler-Silacci et aI]|2006D ? Part of the reason is that 
the brighter disks around early type stars that domi- 
nated high resolution spectral surveys in the past (with 
the ISO-SWS, for example) do not show the same strong 
emission as seen from T Tauri disks. In fact, disks around 
early type (A and B) stars generally lack molecular emis- 
sion strong enough to be detected by the Spitzer IRS 
SH modules. It cannot be ruled out that some Her- 
big Ae/Be stars will have molecular emission, as is in- 
dicated by Figure [M] One apparently exceptional B star 
in this sample shows emissions due to CO2. Further, 
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OH has been detect ed at 3^m in several Herbig stars 
(|Mandell et al.|[2008D . while rov ibrational CO emission 
is common in such s ystems fe.g.JBlak e fc Bo ogerti[200^ 
iBrittain et 311120071 : Ivan der Plas et a l. 2009). The inci- 
dence rate of detectable molecular emission from Herbig 
stars in Spitzer spectra appears to be at least 10 times 
less common than for disks around later type stars, how- 
ever. The reason for this is presently unclear, but could 
be due to a combination of photodestruction of molecules 
by the strong UV fields of the A and B stars, a weakening 
of lines due to physical conditions of the disk, such as a 
reduced gas-to-small-dust ratio relative to T Tauri disks, 
and masking by a strong infrared continuum. 

Further, this lack of molecular emission from Herbig 
Ae/Be stars shows the importance of extending surveys 
with Herschel to spectral types of at least G and F. How- 
ever, the low excitation water lines in the Herschel-PACS 
range are so optically thick that water abundances even 
as low as 10~^° are expected to produce stron g lines 
(|Meiierink et all l2008t [Pontoppidan et all I2OO90 . The 
mid-infrared lines, in contrast, require water abundances 
in the 10"'' — 10~* range, and are therefore important 
tracers of inner disk chemistry. 

The Spitzer data indicate that planetesimals within the 
snow line generally form in a gaseous environment with 
a high water abundance, assuming that vertical trans- 
port is efficient enough to ensure that the disk surface is 
representative of the interior within the midplane snow 
line. This suggests that the formation of oceans on ter- 



restrial planets may not require seeding by a late veneer 
from the outer reaches of a planetary system, such as the 
late heavy bombardm ent event, if water can efficiently 
adsor b to grains (e.g. IDrakd 120051 : IMuralidharan et all 
I2008D . Conversely, the presence of abundant water va- 
por within the snow line strongly hints at a large, unseen 
reservoir of ice beyond ~ 1 AU. 

The mid-IR bands of water and other molecules are 
key tracers of planet formation. Given their ubiquity, 
future observations by the JWST, SOFIA and ground- 
based facilities such as the European Extremely Large 
Telescope (E-ELT), the Thirty Meter Telescope (TMT) 
and the Giant Magellan Telescope (GMT) will have a rich 
list of targets and provide an abundance of constraints 
for chemical models of inner protoplanetary disks. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. Support for this work 
was provided by NASA. Support for KMP was provided 
by NASA through Hubble Fellowship grant #01201.01 
awarded by the Space Telescope Science Institute, which 
is operated by the Association of Universities for Re- 
search in Astronomy, Inc., for NASA, under contract 
NAS 5-26555. Research support for JSC was also pro- 
vided by 6.1 base funding at the Naval Research Labo- 
ratory. The authors would like to acknowledge valuable 
discussions with Uma Gorti, Ilaria Pascucci, and Ewine 
van Dishoeck. 
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